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. We calculate the corrections of orders , and to the Lamb shift of the IS and 

Q \ 2S energy levels of muonic hydrogen (/xp) and muonic deuterium {fid). The nuclear 

structure effects are taken into account in terms of the proton rp and deuteron 
charge radii for the one-photon interaction and by means of the proton and deuteron 
electromagnetic form factors in the case of one- loop amplitudes. The obtained nu- 
merical value of the isotope shift (//d) - {^p) for the splitting (IS" 25) 101003.3495 
^ , meV can be considered as a reliable estimation for corresponding experiment with 

' the accuracy 10-6. The fine structure interval E{IS) - 8E{2S) in muonic hydrogen 

• and muonic deuterium are calculated. 

(N 

£j ■ PACS numbers: 31.30.Jv, 12.20.Ds, 32.10.Fn 

o 

^7^. I. INTRODUCTION 

Oh! 

r-| ! The investigation of the energy spectra of hydrogenic atoms is of great importance for 

high accuracy verification of the Standard Model and derivation more correct values of 
, . ■ fundamental physical constants (the fine structure constant, the masses of the muon and 
^ ! electron, the proton charge radius etc.) 0, Q, ^^^^ years significant interest in 

I this research area is connected with the atom of muonic hydrogen 0, IM IE Q ■ motivated 
by the activity of the experimental study of the {2P — 2S) Lamb shift and hyperfine structure 
in muonic hydrogen. The measurement of the (2P — 25") Lamb shift with the accuracy 30 
ppm in the (fip) would allow us to obtain the value of the proton charge radius with a 
precision 10"^ which is an order of the magnitude better in the comparison with different 
methods including the elastic electron-proton scattering and the value of the (2P — 25") 
Lamb shift in the atom of electronic hydrogen. The measurement of the ground state 
hyperfine splitting in muonic hydrogen with a similar accuracy would allow to determine 
the value of the Zemach radius ^] with a precision 10-3 0,113^ which can be considered as 
a new fundamental parameter regarding to the hydrogen atom. Then it can be used in the 
calculation of new theoretical quantity of the hyperfine splitting in electronic hy drogen and 
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the restriction on the value of the proton polarizability contribution 

There exist another experimental task of the investigation of large fine structure interval 
(IS* — 2S) in muonic hydrogen and the muonic hydrogen - muonic deuterium isotope shift 
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for this splitting 3, 1^ • It would allow to obtain new data about the proton and deuteron 
charge radii. It should be pointed out that both marked intervals are among the most 
precise measured quantities for the atom of electronic hydrogen. Experimental accuracy 
for the measurement of the hydrogen-deuterium isotope shift for the transition (IS* — 25*) 
increased by three orders of the magnitude during last ten years. At present time it is equal 



Aujs = [E{2S) - E{1S)]d - [E{2S) - E{1S)]h = 670 994 334.64(15) kHz. (1) 

The interval (IS* — 25*) for the hydrogen atom was measured with extremely high accuracy 
of order of lO'^ kHz Q: 



Auis-2s{H) = 2 466 061 413 187 103(46) Hz, 6 = 1.8- 10 



^14 



(2) 



Experimental investigations of the energy intervals (1) and (2) in muonic hydrogen are on 
the stage of preliminary preparation. 

Theoretical studies of various contributions to the energy levels of muonic atoms were 
made long ago [l^ |2l| (see also the Refs. in the review article p|). In the last years 
the calculation of different order corrections in the energy spectrum of muonic hydrogen 
was connected for the most part with the (2P — 25") Lamb shift and hyperfine structure 



of S'-states [l^, 0, 0, 2^. In these papers the two-particle interaction operator was 
constructed which gave the contributions of orders and to the interval (2P — 25) and 
hyperfine splitting of the 15* and 25* states. At present the necessity exists for theoretical 
investigation of the corrections of orders a^, and to the 15* and 25* Lamb shift in 
muonic hydrogen and muonic deuterium, the isotope shift (fip) - [fid) for the transition 
(IS* — 2S) and the fine structure interval [E{1S) — 8E{2S)] which are yet unknown. Such 
calculations can initiate the experimental investigations of the fine structure intervals (1) and 
(2) in muonic hydrogen in order to obtain more exact values of such fundamental physical 
constants as the proton and deuteron charge radii and the muon mass. 

In this paper we have derived numerical results for the contributions of orders a^, and 
to the splitting (IS' — 25") and the isotope shift (/xp) - (/irf) for this splitting. In doing 
so, we obtained numerical values for some contributions on the basis of known analytical 
expressions. The most part of the corrections for the energy levels IS* and 25* and the isotope 
shift in muonic hydrogen is determined originally in the integral form which then is used 
for numerical estimations. The dependence on principal quantum number n for the most 
part of the contributions is not trivial that is not reduces to the factor Xjr?. The reason 
has to do with characteristic value of the particle momenta and will be discussed below. 
The aim of this work consists in the calculation of quantum electrodynamic corrections in 
the muonic hydrogen - muonic deuterium isotope shift for the transition (15* — 25) and the 
fine structure interval \E{\S^ — 8E{2S^\ both for (/ip) and (/id) with the accuracy 10~^ 
which can be considered as a proper guide for the future experiments and the retrieval of 
the information about the proton and deuteron charge radii and the muon mass. 

Fine structure of the energy levels of hydrogenic atoms has been studied for a long time 
on the basis of different approaches [l|, i26., 23\ ■ For the case of S'-states of hydrogen-hke 
atoms consisting of the particles with masses m\ and it can be written with the accuracy 
0((Za)^) as follows: 



E„ = mi -t- m2 



2n2 



2n3 



+ 



An Amim2n 



(3) 
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fip{lS) : 1 043 927 826 470.3586 meV; ^23) : 1 043 929 722 866.0601 meV 
fxd{lS) : 1 981 268 455 762.7537 meV; fxd{2S) : 1 981 270 453 188.8081 meV 



The error of relative order in Eq.(3) which is determined by the uncertainties of the fine 
structure constant a and the particle masses reaches the value 10~^. Nevertheless, we 
represented the numerical values in Eq.(3) with the precision 0.0001 meV because it is 
important when we consider various fine structure intervals in the energy spectrum. The 
following values of fundamental physical constants are used: a^^ = 137.03599976(50), 
= 0.105658357 (5) GeV, = 0.938271998 (38) GeV, = 1.875612762 (75) GeV 0. The 
contribution of Eq.(3) to the (/xp) - (fid) isotope shift for the splitting (IS' — 25") is the crucial 
(see the Table I). In addition, there exist a number of important effects both electromagnetic 
and strong interactions which must be taken into account for the determination of the isotope 
shift value with the accuracy up to terms of order a^. 



II. EFFECTS OF ONE-LOOP AND TWO-LOOP VACUUM POLARIZATION IN 

ONE-PHOTON INTERACTION 



Our calculations of various energy levels of hydrogen-like atoms are performed within 
the framework of the quasipotential approach in which a bound state of two particles is 
described by the Shroedinger-type equation 2^ 2^: 



M 



2hr 2/i 



^Af(p) 



(iq 

(2^ 



■r(p,q,M)V^M(q) 



(4) 



where 



El 



m: 



'21 



Hji = E1E2/M is the relativistic reduced mass, M = Ei + E2 is the mass of the bound state. 
The quasipotential in Eq.(4) is constructed in quantum electrodynamics by the perturbation 
theory with the use of the two-particle scattering amplitude T projected onto the positive 
frequency states outside the mass surface at zero relative energies of the particles: 



(5) 
(6) 



The initial approximation of the quasipotential V{p, q, M) for a bound system was se- 
lected in the form of the usual Coulomb potential. The increase in the lepton mass in muonic 
hydrogen compared with its electronic counterpart decreases the radius of the Bohr orbit in 
the (/ip). As a result, the Compton wave length of the electron and the radius of the Bohr 
orbit become commensurable: 



h 



0.737384, 



(7) 



where me is the mass of the electron, is the reduced mass in the (/xp) atom. This sub- 
stantially enhances the role played by vacuum polarization effects in the energy spectrum 
of muonic hydrogen [3^. Corrections of the one-loop and two-loop vacuum polarization to 
the quasipotential of one-photon interaction are shown in Fig.l. 
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abed 
FIG. 1: Effects of one-loop and two-loop vacuum polarization in one-photon interaction. 

To determine the contribution of diagram (a) in Fig.l (electronic vacuum polarization) to 
the particle interaction operator we must perform the following substitution in the photon 
propagator js^: 

\_ a p 7^^(2.^ + 1) 
k^^SirJi s4(P + 4m2s2) ■ 

If 

(electronic hydrogen, is the reduced mass of two particles in the hydrogen atom), then, 
ignoring the first term in the denominator in the right-hand side of Eq. (8) we obtain: 

— a/lSvrmg. 

However, if 

(muonic hydrogen, mi is the mass of the muon), then fia and rrif. are values of one order 
and we cannot expand the denominator in Eq.(8) in a. With muonic hydrogen, we must 
construct the particle interaction operator in the one-photon approximation using exact 
expression (8). Further, we take into account that the electron vacuum polarization gives 
the contributions of orders a^, and to the energy spectrum of ^-states in the (/xp) 
atom. 

The modification of the Coulomb potential 

V^{k) = -Ze^k^ 

caused by the vacuum polarization (VP) is determined taking into account (8) by the fol- 
lowing expression in the momentum representation (30l |: 



The Fourier transform of Eq.(9) gives the corresponding operator in the coordinate repre- 
sentation: 
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FIG. 2: The Wichmann-Kroll correction. The wave hne denotes the Coulomb photon. 



The potential (10) allows us to obtain the correction of the electron vacuum polarization of 
order to the energy levels of IS* and 25'-states in muonic hydrogen. Accounting exact 
expressions of the wave functions for 15* and 25'-states 

V^ioo(r) = , %l)2mir) = -^=6 M 1 - — , W = fiZa, (11) 



we represent this correction in the form: 

/i(Za)2a foo 1 rue^ ^/^t^(2^2 _^ 1) 



(12) 

AE,^yp{2S) = - ^^ I piOd^ ( 3^ - 377V + T^ZJTTV I ' ^2(0 = 1 + (13) 



The electron vacuum polarization effects are very sensitive to the bound state structure 
because the characteristic momentum of the particles in muonic hydrogen atom is {fiZa). 
So, the contribution of the vacuum polarization amplitudes in one-photon interaction cannot 
be expressed by a simple factor |'0'-^(O)p ~ but it is dependent on principal quantum 
number n in more complicated form which is represented by expressions (12) and (13). 
Numerical values of the electron vacuum polarization contribution to the 15* and 25* Lamb 
shift in muonic hydrogen and muonic deuterium differ only due to reduced mass of two 
particles: 

_ f fip{lS) : - 1898.8379 meV; pp{2S) : - 219.5849 meV 
^'^'^^ ~ \lid{lS) : - 2129.2820 me\/; /ici(25) : - 245.3205 me\/ ^ ' 

The contribution of the muon vacuum polarization (MVP) can be found also with the use of 
Eq.(9) in which the replacement irte — > mi must be carried out. This correction of order 
to the energy spectrum is included in Table I together with the muon self energy correction 
(MSE). 

The diagram in Fig. 2 has the form of one-loop vacuum polarization. It determines the 
correction to the energy spectrum of order known as the Wichmann-Kroll correction 
(31 . 3^ . The corresponding potential was represented in the integral form: 



vrr Jo C 1^ Jq 



2 I K 

I- - ■ • /•= / - - WK, 



X 



(15) 
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Exact expression for the spectral function is written in Ref . P. Isli 3^ . Averaging the 
potential (15) over the Coulomb wave functions (11) we expressed the corrections to the 
S'-states as follows: 



IT 



dC 



71 

12 



CV.iO 



C - ie{c - 1) + / dx^e - x'f 



2fWK, 



X 



dC 



+ 



X 



71 

"12 



2vr Jo \pliO pHO 2pt(C), 

c 



(16) 
(17) 



2fWK, 



X 



Numerical values of Eqs. (16) and (17) are included in Table I. 

Let us consider the modification of the Coulomb potential by virtue of two-loop vacuum 
polarization (see Fig.l(b,c,d)). The contribution of the diagram (b) in Fig.l containing two 
sequential electronic loops can be obtained by applying the replacement (8) two times in 
the photon propagator. In the coordinate representation the required interaction operator 
of the particles has the form: 



^17, vp-vpi''^) 



a 

9^ 



/OO /•OO / 



Za 



r J {i"^ — 1]"^ 



It gives the following results for the energy spectrum: 

jj,a'^{Zay 



AE 



I7, VP-VP 



(IS) 



97r2 



e 



(^2 _ ^2) \pl[^) pl[ri) 



(18) 

(19) 



X 



vp-vpi'^S) 
1 2 3 



j/ip: — 1.8816 mey 
"j/irf: - 2.1871 mel/ 

l2a'^{Zay 



187r2 
/ 1 



piOd^ / p{v)dv 



+ 



X 



(20) 



fxp: -0.2A26meV 
fxd: - 0.2811 mel/' 



The contributions of the diagrams (b,c) in Fig.l which are determined by the second 
order polarization operator can be calculated after the following replacement in the photon 
propagator [33i] : 



1 



a 



71 



2 rl 



f{v) 



Lin 



16^ ' 4 



Ami + ^^(1 ^ 



dv={-] - 



dv 



7c J 3 Jo + fc2(l — t>2] 



X 



(21) 



In 



+ 



+ 2Li2 
3 



1-v 
1 + v 



3, 1 + v 

H — m m 

2 1-v 2 



, 1 + v ■ 
m m V 



+ 



-v(3 — f ^) In — 

2 ^ ' 4 



2v{'i - v'^)\\iv 



+ -17(5-3^2 
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The contribution value can then conveniently be calculated in the coordinate representation 
with the reduction of the interparticle interaction potential to the form: 



_ _2_Za /«y ^1 f{v)dv 



3 r 



(1 - v"") 



(22) 



The operator (22) gives the following corrections to the Lamb shift of ^'-levels in muonic 
hydrogen and deuterium: 



AE 



I7, 2-loop VP 



2 ^2W"^ 

jia [Za) 



37r2 



ml Jo 



f{v)dv 



V me J 



HP : - 12.6144 meV 
lid : - 14.0141 meV 



AE 



I7, 2-loop VP 



{2S) 



127r2 



1 W"^ 



f{v)dv 



(1 + 



2me 



Pz{v) 2pl{v) 



(23) 



(24) 



HP 



1.4112 mey 



//d: - 1.5606 mey 



1 



1 + 



2me 



Note that, as we determine contributions to the energy spectrum numerically, the corre- 
sponding results are given with an accuracy of 0.0001 meV. 



III. THREE-LOOP VACUUM POLARIZATION IN ONE-PHOTON 

INTERACTION 

In the fifth order over a there is the contribution of three-loop amplitudes of the vacuum 
polarization in the one-photon interaction (see diagrams (a)-(b) in Fig. 3). The contribution 
of the diagram (a) in Fig.3 to the potential takes the form: 



VP-VP-Vp{''^) 



Za 



r (37r) 



X 



-2meC,i' _ 



+ e 



-y'J poo rco rco 



X 



(25) 



+ e 



ie-ew-c) ic-ew-e) ie-v'w-v') 

The contribution of the diagram (b) in Fig.3 can be written also in the integral form 
using expressions (8) and (21): 



c 



AiJia^{Za) f°° 



VP-2-loop VP 



97r3 



rj r 



-2mer]r V 

rf — 



^~2me^r _ 



(26) 

The corrections to the energy spectrum of the {/ip) and {/id) atoms corresponding to the 
interaction operators (25) and (26) are the following: 



AE^ 



VP-VP-VP 



{IS) 



lia^{Za) 
277r3 



X 



+ 



; piOd^l P{v)dvl piOdCx 



(27) 



+ 



8 



G 



G 



FIG. 3: Effects of three-loop vacuum polarization in one-photon interaction (a,b) and in the third 
order of the perturbation theory (c). 



VP-VP-VP 



{2S) 



j fxp: - 0.0029 meV 
[fid: - 0.0034 me 



X • 



547r3 
1 



P{v)dvl p(CKx (28) 



VP-2-loop VP 




0.0003 meV 
0.0004 meV 

e 



+ 



pm 



AE- 



VP-2-loop VP 



(25) 



_ j fxp: - 0.0223 meV 
~ [fid: - 0.0251 meV 



fiv) dr] 



(29) 



pHv) pUv) 



(30) 



HP : - 0.0030 meV 
fid: - 0.0035 mel/' 



[pm pm 

There is need to make the substitution v = y/rj'^ — l/rj in the function f{v) which is defined 
by Eq.(21). There exist also several diagrams which represent three-loop corrections to the 
polarization operator. They were calculated originally for the Lamb shift (2P — 25") in the 
papers of the Kinoshita and Nio 3^. The greatest contribution to the energy levels is 
determined by the diagrams of the sixth order vacuum polarization with the one electron 
loop (the contribution 11'^^^) |13)- The estimation of their value to the 15* and 25* Lamb 
shift is presented in Table L 



9 



IV. THE VACUUM POLARIZATION AND RELATIVISTIC EFFECTS 



To calculate the energy spectra of the ^'-states of muonic hydrogen with a precision up to 
terms of order we have to construct the quasipotential on the basis of relations (5) and 
(6) which catch correctly the relativistic effects of necessary order (the Breit Hamiltonian 
AVb)- Accounting the electron vacuum polarization effects Pachucki has derived the Breit 
Hamiltonian AV^^ in Ref. j2l|. Omitting the spin-dependent terms of the interaction 
operator let us write the two-body Breit Hamiltonians AVg and /SV^^: 



AV, 



B 



-P 



1 1 
+ — 



TxZa 



ml 



1 1 



a 
Stt Ji 



Za 



ml 



1 1 

+ 



5(r)- 



Za 



2mim2r 



2 I nrjp.pj 



2 

mf 



ml 



Zaml(^ -2mfr/i ^ \ Za e 



7r5(r) — 



(31) 
(32) 



5,, + ^(l + 2meer) 



Pr 



mim2r 2mim2 r 

In the first order of the perturbation theory (FOPT) the Hamiltonian (32) gives the following 
contribution after averaging over the Coulomb wave functions (11): 

. j^reiyp ^ , I . ,.yP| , ^ / M^S) : 0.1962 meV; iip{2S) : 0.0249 meV 
AE,^^ -< te^lAV^B >- . Q 2515 meV; ^^d{2S) : 0.0322 meV 

(33) 

These corrections are of order a{ZaY. Second-order perturbation theory (SOFT) cor- 
rections to the energy spectrum are determined by the reduced Coulomb Green function 
(RCGF) '35'], whose partial expansion is written as: 

G„(r,r') =Y.gni{ry)YU^)Y:^{r^). (34) 

The radial function gni{,i^) 1^') was obtained in Ref. 3^| in the form of the Sturm expansion in 
the Laguerre polynomials. For the IS" and 2S'-states these functions are the following: 



9io[r,r 



-2fi^Za 



^ \t2 m(m-l) 2 2 2) 



(35) 



m{m — 2) \2 dx dx\ 



L\{x)Ll{x') 



m=l,mf^2 

where x = fiZar, are the usual Laguerre polynomials defined as: 



e 2- 



w. 



(36) 



(37) 



As a several quasipotential terms contain 5(r), we must know Gn{r,0). The correspond- 
ing expression for the reduced Coulomb Green function was obtained in Ref. [37l| using the 
Hoestler representation for the Coulomb Green function and subtracting the pole term. This 

gave 



Zafi^ 2e-^'/2 



X 



2x{\nx + C) + - 5x - 2 



(38) 
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ZoLiP- e~^/^ r 1 
G25(r, 0) = 4x(x - 2)(lnx + C) + - ISx^ + 6x + 4 , (39) 

where C = 0.5772... is the Euler constant. One-loop vacuum polarization gives the contri- 
bution in the second order of the perturbation theory which is specified by the relation: 



SOFT 



2 E 



EG -E£ 



(40) 



The matrix element of the operator is expressed by means of the substitution p^/4//^= 
{Hq + Za/r){HQ + Za/r) {Hq — p^/2// - Za/r). After that we employ the sequence of 
algebraic transformations: 



<^n|— ^0 2^ FC_ PC ^VPWn >- 



(41) 



2^ P(7 rr7 Vpl^n -> ■ 



Integrating the RCGF (34) over the coordinates r and r' we obtain the sums which are 
calculated accounting the range of a variable C (1 -^ oo) and 1/(1 -|- W/rrieC) < 1- Typical 
matrix element for the IS* correction is the following: 



5 



X 

2 

+ 



X 

2" 



+ E 

lnpi(0 



-^m-l(^)-^m-l(^0 

m(m — 1) 



= (42) 



4p?(0 2pf(0 p?(0 ' 

[pi(o-ir-^ _ 1 1 



E 



lnpi(0 



Omitting other numerous intermediate analytical expressions we write the summary numer- 
ical value of the second order of the perturbation theory contribution (40) for the l^S" and 
2S'-levels: 



^^reiyp _ \ ^^P{'^S) : - 0.2644 meV] fip{2S) : - 0.0559 meV 
^' ~ \fid{lS) : - 0.3194 meV; i^d{2S) : - 0.0696 meV 



(43) 



V. TWO-LOOP AND THREE-LOOP VACUUM POLARIZATION IN THE 
SECOND ORDER OF THE PERTURBATION THEORY 



Two-loop vacuum polarization gives the correction in the second order of the perturbation 
theory (see the diagram (a) in Fig.4): 



c 



(44) 



The calculation of this matrix element can be done with the relations (10), (35) and (36). 
As a result we obtain the corrections of order a^{ZaY for the 15" and 25" states: 



lJia^{Za)" 



2pl{C)pl{r,) 



(45) 
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FIG. 4: Corrections of two-loop and three-loop vacuum polarization in the second order of the 
perturbation theory. 



1 



1 



+ 



[pi(0-i][pi(^)-i] 



1 



In 



[pi(0+pi(^)-i] 

9 6 3 



(46) 



+ 



9 



+ 



153 



39 



MiOpUv) pUOpUv) pHOpUv) '2pmpi{v) SpiiOpiiv) ^pUOpiiv) 

33 6 39 17 6 3 



+ 



+ 



+ 



4:pl{0pi{v) pUOpUv) '^pUOpUv) pUOpUv) pHOpUv) pUOpHv) 

3 z\6 - 5z) 



+ 



33 



6 



"^pUOpUv) pUOpHv) 

^"^(6 — 5z) 



+ 



z\4 - 3z) 



z\4 - 3z) 



-30 + 49^-21^2) 



2ai(0ai(r^)(^ - iyp,{r^) 2ai(0ai(r/)(z - 1)^ 2ai(0ai(r^)(^ - 1)^ 44(^4^^ - If 



-18 + 27^-11^2 



^^(-18 + 27^-11^^ 



^^(-10 + 13^-5^2) , 
+ w + ln(l - ^) X 



4ai(0ai(r/)(^-l)3p2(0 4ai(C)ai(r;)(;. - 1)2^2(77) 4ai(0ai(r^)(^ - 1) 



X 



2^3 



3z^ 



3z' 



9z^ 



+ 



3z^ 



+ 



3z^ 



2amaUv) 2amaliv)P2iO 2aliOaliv)P2iv) 2ai(0ai(r;) J ' 

where 02(0 = P2iO - 1= ^ = a2(0«2(?7)/p2(C)P2(?7)- 

Three-loop vacuum polarization contributions in the second order of the perturbation 
theory shown in Fig.4(b,c,d) are of order a^{Zay. For their calculation we must use again 
the analytical expression for the one-loop and two-loop polarization operator. Fulfilling the 
integration over the wave function coordinates we represent the necessary corrections in the 
form: 



A jr^vp-vpyp 



277r3 



p{r])dr] 



n = ^-J- ^ + [pii^)-M[pi{C)-i] 

^^^'''^ 2 pi(0 pi(C) Pi(0+Pi(C)-i 



1 {e-v')J^ pKC) ' pKO 

piiO+PiiO - 1 



Pliv) 



Pi(0-Pi(C) 



(47) 
(48) 
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AE 



vp-vpyp 

SOFT 



{2S) 



277r3 



/oo poo 
Piv)dvl PiOdC 



[efi^,c)-v'fivx)] 



where the function f is defined earher by the Eq.(46) 



A j-^2-ioop vpyp 

^-^SOPT 



(IS) 



1 



1 



pi(0?i(^) 



9vr2 Jo ^^''^'^l2p?(Ogf(tO Pimiv 



+ 



In 



AE 



PiiO + Qiiv) + ^ PiiO + Qiiv) 

Afia^Zaf /•! f{v)dv 



qi[v 



1 + 



2-ioop vpyp 

SOPT 



(25) 



97r3 



1 



Ji 



where P2{v) = 1 + 2me/Wy/l — in the function f{^,v), a2{v) = 2me/W^/l 
energy terms (45)-(51) give the contribution to the (fip) - (fid) isotope shift: 

ypyp 



AE 



SOPT,IS 



0.3125 meV, 



(49) 
(50) 



(51) 
The 

(52) 



and the individual values for the S-states are depicted in Table I. The analysis of the three- 
loop vacuum polarization contribution in the third order of the perturbation theory (see 
the diagram (c) in Fig. 3) shows that it is one order of the magnitude smaller than the 
contribution of the amplitude shown in Fig.4(c) and can be neglected. 



VI. MUON SELF-ENERGY, MUON VACUUM POLARIZATION AND RECOIL 

CORRECTIONS 

To estimate the corrections of the muon self-energy (MSE) and muon vacuum polarization 
(MVP) of order a{Za)^ we used well-known analytical results for the hydrogen atom. In 
the case of the IS' and 2S'-states these contributions are the following jH: 



AE 



MVPMSE 



AE 



MVP,MSE 



{IS) 



{2S) 



IT 



2 



rrti 



3 j2{Zay 



57r m 



4 J mi 
3 fi{Zay 



-lnA;o(l,0) + — 
3 ' ^ 45 



--lnA;o(2,0) + — 
3 ' ^ 45 



fip : 5.1180 meV 

fid : 5.9395 meV 
(53) 

fip : 0.6543 meV 
fid: 0.7594 mel/' 
(54) 



where lnA;o(^;0 is the Bethe logarithm: 

In ko{l, 0) = 2.984128555765498, In ko{2, 0) = 2.811769893120563. 



(55) 



The correction determined as a sum of radiative insertions in the lepton line is known also 
in analytical form. We represent its value despite the fact that the order of this correction 
is a{Za)^ because the corresponding coefficient is sufficiently large (of order 10) [l|: 



nil \ 96 



3mT 



96 



f fip : 0.0355 meV 
[fid: OMUmeV 

fip : 0.0044 meV 
fid: 0.0052 meV 



(56) 
(57) 
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FIG. 5: Radiative corrections accounting electron vacuum polarization. 



Following the Pachucki paper |21|, let us estimate also the contributions of two diagrams of 
the sixth order over a enhanced by the In a. The diagram (a) in Fig. 5 of the muon radiative 
correction with electron vacuum polarization gives the following contributions: 



rad+VP 



(IS) 



2 1 

+ 



Piio pm pm 



(58) 



X 



16 9 3 



fip : 0.0061 meV 
jid : 0.0073 meV 



9 



In^/ piOd^x 



(59) 



+ 



■ ftp : 0.0010 meV 
l^d : 0.0012 meV 



.3P2(o 2pi(o pm ^piio pm. 

The diagram (b) in Fig. 5 gives the contribution to the energy spectrum which can be 
expressed in terms of the slope of the Dirac form factor F[ and the Pauli form factor F2: 



AE. 



2,rad+VP 



nS) 



rrii 



4mfF{(0) + F2(0) 



(60) 



where the form factors F[ and F2 were calculated in Ref. 

2 



miFiiO) 



1 , n mi 29 , mi 1 ^.^^ 395 

- In^ — In — + -C(2) + 

9 me 108 me 9^ ^ 1296 



^2(0) 



2rl^^mi 25' 



rrip 



36 



Then the correction to the Lamb shift of the S-levels takes the form: 

, ^, fi^a^Zay/A.mi 20, mi 4^,, 85 
AE2,rad+vp{nS) = ^\ J In^ ^ - - In ^ + -( 2 + — 
mfn'^TT^ V9 me 27 me 9 162 



(61) 
(62) 

(63) 
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In this section we included yet known analytical recoil correction of order iZaf' which is 
determined by two-photon exchange diagrams js^: 



where 



(Za)V 
TTn^mim2 

2 

2 2 
1712 " 'f^l 



n 



2 8 17 

-- In(Za) - 2 In koin, 0) - g - 



mi 2 1 "^2 



In In — 



1 + ^ + ... + - 

2 n 



1 

2n 



Numerical values of the corrections (63) and (64) are presented in Table I. 



(64) 



(65) 



VII. EFFECTS OF NUCLEAR STRUCTURE, POLARIZABILITY AND 

VACUUM POLARIZATION 

In the energy spectrum of muonic hydrogen the important role belongs to strong interac- 
tions which are connected with the distributions of electric charge and magnetic moment of 
the nucleus. Namely these corrections together with the mass differences of the isotopes lead 
to the isotope shift in the system (/ip) - (/id). In the leading order {ZaY the nuclear struc- 
ture effects are determined by the differential parameter of the electric charge distribution 
known as the nucleus charge radius r^v- To calculate one- loop corrections where the nucleus 
structure effects are essential we must know the form of the nucleus electromagnetic form 
factors. The contribution of the nuclear structure effects was studied both for the hyperfine 
structure and the Lamb shift in hydrogenic atoms in Refs.[l|, Hll, 113, 113]. The leading 
order nuclear structure correction of order {ZaY to the Lamb shift of the ^'-levels in muonic 
hydrogen has the form (see the diagram (a) in Fig.6): 

where is the mean-square proton charge radius. Numerical values of the correction (66) 
for the levels with n = 1 and = 2 (at rp=0.891 fm ^) are presented in Table I. They give 
essential relative order contribution, so, the decrease of the error in the value of the proton 
charge radius is extremely important task if we desire to obtain more precise theoretical 
values of the Lamb shift of the S'-states, different fine structure intervals and the isotope 
shifts. In the case of the muonic deuterium the deuteron charge radius rrf=2.094 fm is used 

Two-photon exchange amplitudes shown in Fig. 7 give the nuclear structure corrections 
of order {Za)^. In this case the quasipotential of two-photon interaction can be found by 
means of the relations (5) and (6) i^]. Corresponding correction to the energy spectrum 
takes the form of one- dimensional integral: 

^i\Zaf dk 



AE 



nS)^-^^ri^V(k), (67) 

nn^ Jo k 



str,{Za)5 

V(t) ^ - ^> + 8^.|ff ) ^ , (68) 

mim2 m2\mi + m2)k 2m'lm2 
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FIG. 6: The nuclear structure and vacuum polarization corrections. 



X 



2{F^ - l){mi + ml) + 4FiF2m^ + 3F|m? 



+ 



2m\m2{m'i — m^A; 
lQm\ml{Fl - I] 



+ 4m2mi 
2m2(mf — m^)/;; 



2(Fi2-l)+4FiF2 + 3F2' 



8mlFiF2 + 



J 

16m^(F2 - 1^ 



To carry out numerical integration in Eq.(67) we used the parameterization of the Dirac Fi 
and Pauh F2 proton form factors obtained in Ref . Jill] . For the muonic deuterium the similar 
contributions are obtained in Ref. |25|. The fifth order over a contribution connected with 
the nucleus structure is given also by the effects of the electron vacuum polarization which 
are shown in Fig.6(b,c). The particle interaction operator corresponding the amplitude in 
Fig. 6(b) is determined by the following expression: 



stryp[r ) 



2a{Za)r 



N 



9 



(69) 



The matrix elements of the operator (69) over the functions (11) lead to numerical data: 

2 



AEstrypilS) = -a(Z«)Vr^ / 

AE,tryp{2S) = 4a(^«)Vr^ / 
06 Ji 



Pirn 



w 

4me^ 
W ' 



fip : 0.1991 meV 
fjid : 1.4155 meV^ 



(70) 



2c2 



4mf^ 



+ 



W \ P2{0 2pi(0 



(71) 



_ j pp: 0.0257 meV 
~ : 0.1824 meV 

There exist also the contribution of the electron vacuum polarization and nuclear structure 
in the second order of the perturbation theory (see the diagram (c) in Fig. 6) which is 
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FIG. 7: The nuclear structure corrections of order (Za)^. Bold circle on the diagram denotes the 
nucleus vertex operator. 



determined by the reduced Coulomb Green function Gn{r, 0) (38) and (39). In this case the 
contributions to the shifts of the S'-levels are: 



2 

97r^ 



2-3pi(0-2p?(0-2pi(0 lnpi(0 



AE. 



str,VP;SOPT 



{2S) 



36n 



HP : 0.1242 meV 
fid : 0.8913 meV 

Pirn 



(72) 



-12 + 23p2(0-8p^(0-4p^(0 + 4p^(0 + 

(73) 



' fip : 0.0126 meV 
fid: 0.0898 mel/' 



The strong interaction contribution to the energy spectrum of the (fip) and (fid) comes 
not only from the nucleus structure but also from the nucleus polarizability. For muonic deu- 
terium atom we used analytical expression of the deuteron polarizability correction obtained 
in Ref.Hi: 



AEd,poiinS) 



a{Zay 



-fi mi 



1 



mi 20 



5«, In — + —-/?, In 1.24 



8/ 



mi 



(74) 



where the deuteron bound energy I = K^/m2, and k = 45.7 MeV, ad = 0.635 fm^ and Pd = 
0.073 fm^ are the electric and magnetic deuteron polarizabilities. For muonic hydrogen atom 
the proton polarizability contribution to the Lamb shift of the S'-states can be calculated 
employing the expression |i22|, |23i] : 



AEp^poi{nS) 



16/i^(Za)^mi dk 



k Jo 



dy 



sm 



1^0 {k'^ + 4mf cos^ 0) (y^ + k'^ cos^ 



X 



^1 + k^ /y ) cos' 



3^,5 



(75) 



X [(1 + 2 cos^ 0) ' + sin' 0]i^2(l/, k^) + 

1 + R[y,k^) 



7rn-^mim2 Jo 



h{k^)(3{e)kdk, 
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where F2 is the structure function of the lepton-nucleon scattering, R = a^/ax is the ratio 
of the cross sections for the absorption of longitudinally and transversely polarized photons 
by hadrons, uq is the threshold for the vr-meson production. The subtraction term in Eq.(75) 
includes the function 

where j3p = 2.1(0.9) x 10~^ fm^ is the proton magnetic polarizability. The contribution 
of internal deuteron polarizability can be calculated also on the basis of Eq.(75) ji^. It is 
included in total value together with the correction (74) which is dominating. 

Hadronic vacuum polarization contribution of order a{ZaY to the shift of S-levels in 
muonic hydrogen which was studied in Ref. Bill can be presented in the form: 

AEHVp{nS) = / , (77) 

where the spectral function p^{s) is expressed through the cross section of e^e~ annihilation 
into hadrons: 

, a^{e+e- hadrons) 

P ^'^-3.a,,(e+e--./.V-)' ^ ^ 

Dividing the whole integration region over s into the intervals where the cross section of 
e"''e~-annihilation into hadrons is known from the experiment ^] we can carry out numerical 
integration in Eq.(77). Total numerical results are included in Table I. 



VIII. CONCLUSION 

In this paper we calculate different quantum electrodynamic corrections, effects of the 
nucleus structure and polarizability, hadronic vacuum polarization to the Lamb shift of 
the IS* and 25* energy levels in muonic hydrogen, muonic deuterium and isotope shift (/xp) 

- (yud) for the splitting (IS* — 25*). The main goal of the work is to summarize various 
corrections to the fine structure interval (15 — 25') in the (/xp), (/i(i) and the muonic hydrogen 

- muonic deuterium isotope shift for this splitting which would allow to obtain reliable 
theoretical values for the indicated quantities. We consider the corrections of orders a^, 
a^, and also several important contributions of order enhanced by the In a. In our 
calculation we took into account that the ratio [lajnif. is close to one and centred special 
attention on the effects of the electron vacuum polarization. Numerical values of obtained 
contributions are presented in Table I. We included also in Table I the references on many 
papers where analytical or numerical calculations of some corrections are performed despite 
the fact that particular numerical results for the 15 and 25-levels are absent in these papers. 
For the comparison of obtained results with the calculations carried out by other authors 
we used commonly the review article ^ which accumulate recent advances in the physics of 
the energy spectra of simple atomic systems and contain detailed references on the earlier 
performed investigations. 

Let us point out some peculiarities of performed calculations: 

1. The effects of the vacuum polarization for the considered muonic atoms have played 
significant role. They lead to the modification of the Breit two-particle interaction operator 
and give the corrections in the energy spectra up to fifth order in a. 
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2. The nuclear structure effects in tlie energy spectrum of tlie 5'-states are expressed in 
terms of the nucleus (the proton and deuteron) charge radius in the leading order {Za)'^ 
and in the next to leading order (Za)^ for the one-loop amplitudes by means of the nucleus 
electromagnetic form factors. 

3. The estimation of the nuclear polarizability contributions is performed on the basis of 
Eq.(74) obtained by the Mil'shtein, Petrosyan and Khriplovich and Eq.(75) in Refs. jl^, . 
The nuclear structure and polarizability contributions lead to largest theoretical uncertainty 
for the transition (15* — 2S) and the isotope shift (fip) - (fid). 

Total numerical values for the energy levels 15* and 25* in muonic hydrogen and muonic 
deuterium presented in Table I and also the values of the fine structure interval (15* — 2S), 
isotope shift [fip) - {^d) can be considered as a proper estimation for the future experiments 
with these muonic atoms. The values of the corrections are obtained with the accuracy 
0.0001 meV. For the splitting (15* — 25*) the error of theoretical result is determined by a 
number of the factors. The uncertainties of fundamental parameters (fine structure constant, 
the proton, deuteron and muon masses) amount the value of order 10~^. The quantum 
electro dynamic corrections of higher order give theoretical error of order 10~®. The 
largest contribution to the theoretical error is connected with the uncertainties of the proton 
and deuteron charge radii. Their relative contributions reaches the value of order 10~^ (we 
used the values of the proton and deuteron charge radii: Vp = 0.891 fm, = 2.094 fm). 
Further improvements of theoretical results presented in Table I are related primarily with 
the nuclear structure and polarizability corrections. So, for the comparison of performed 
calculations with future experimental data a major interest is connected with the quantity 
of fine structure interval [£'(15') — 8£'(2S')] which has not contain the nuclear structure and 
polarizability corrections of the leading order (Za)^. Numerical values for this fine structure 
interval in muonic hydrogen and muonic deuterium are the following: 

fip : E{1S) - 8E{2S) = -7 307 509 956 602.3099 meV, (79) 

fid : E{1S) - 8E{2S) = -13 868 895 169 916.6917 meV. 

The relative value of that part of theoretical error in Eq.(79) which is determined by the 
QED corrections of higher order is extremely small (of order 10~^^). 

The muonic hydrogen - muonic deuterium isotope shift for the splitting (IS* — 25") is 
considered to be among the most important characteristics of the energy spectra of muonic 
atoms. The differences in wavelengths of lines emitted by isotopes of the same element 
can arise either as a result of the differences in the masses of the isotopes or on account 
of differences in the nuclear charge distributions. Whereas the masses of the proton and 
deuteron are determined at present with sufficiently high accuracy, the nuclear structure 
parameters are known less precisely. Using the results of fulfilled calculations we can express 
the difference of the deuteron and proton charge radii in terms of isotope shift: 



rl 12 



AEf^ - AETs') , (80) 



where theoretical value AEj^ doesn't contain the nuclear structure correction of order {Za)^. 
So, the measurement of the isotope shift (np) — (fid) will allow on the one hand to perform 
additional test of the QED and on the other hand to obtain more accurate value of the 
deuteron charge radius from the relation (80) after the extraction the proton charge radius 
in the experiment at Paul Scherrer Institute (PSI) 01 ■ 
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TABLE I: Corrections of order a^, and ct' to the 15" and 25" Lamb shifts in muonic hydrogen 
and muonic deuterium and isotope shift ^Eis- 



The contribution 
to the energy of the atom 


fip, meV 


fid, meV 


AEis 


Ret. 


IS 


2S 


IS 


2S 


1 


2 


3 


4 


5 


6 


7 


Fine structure 
formula 

En = 1711+1712- 

txiZa)^ L 3 1 n'' ] 
2rv^ in ' 4nmim2 


1043927826 
~470.3586 


1043929722 
~866.0601 


1981268455 
~762.7537 


1981270453 
~188.8081 


101030.3530 


(3) 

UJ 


One-loop VP 
in I7 interaction 
of order a{Zaf' 


-1898.8379 


-219.5849 


-2129.2820 


-245.3205 


204.7085 


(12) 
(13) 
(14) 


Wichmann-Kroll 
correction of 
order a{Za)^ 


0.0114 


0.0012 


0.0126 


0.0014 


-0.0010 


(16) 

u 


Two-loop VP (VP- VP) 
in I7 interaction 
of order a^(Za)^ 


-1.8816 


-0.2426 


-2.1871 


-0.2811 


0.2616 


(19) 
(20) 


Two-loop VP (2-loop) 
in I7 interaction 
of order a'^{Za)'^ 


-12.6144 


-1.4112 


-14.0141 


-1.5606 


1.2476 


(23) 
(24) 


Three-loop VP (VP- 
VP- VP) in I7 interaction 
of order a^(Za)^ 


-0.0029 


-0.0003 


-0.0034 


-0.0004 


0.0004 


[34] 
(27) 
(28) 


Three-loop VP (VP- 
-2-loop) in I7 interaction 
of order a^{Za)^ 


-0.0223 


-0.0030 


-0.0251 


-0.0035 


0.0023 


m 

(29) 
(30) 


Three-loop VP (H^^^)) 
in I7 interaction 
of order a^{Za)'^ 


-0.0340 


-0.0045 


-0.0380 


-0.0050 


0.0035 


[34] 


Relativistic and 
VP correction in 
FOPT of order a^{Zaf 


0.1962 


0.0249 


0.2515 


0.0322 


-0.0480 


(33) 


Relativistic and VP 
effects in SOPT 
of order a^(Za)^ 


-0.2644 


-0.0559 


-0.3194 


-0.0696 


0.0413 


(43) 


Two-loop VP in 
SOPT of order a'^{Zaf 


-2.0343 


-0.1532 


-2.3675 


-0.1750 


0.3114 


(45) 
(46) 


Three-loop VP 
(VP- VP, VP) in SOPT 
of order a^{Za)^ 


-0.0061 


-0.0002 


-0.0073 


-0.0005 


0.0009 


(47) 
(49) 
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Table I (continued). 



1 


2 


3 


4 


5 


6 


7 


Three-loop VP 
(2-loop VP, VP) in 
SOPT of order a^{Zaf 


-0.0059 


-0.0016 


-0.0069 


-0.0021 


0.0005 


(50) 
(51) 


Muon SE and VP 
of order a{Za)'^ 


5.1180 


0.6543 


5.9395 


0.7594 


-0.7164 


(53) 
(54) 


Radiative 
corrections 
of order a{Za)^ 


0.0355 


0.0044 


0.0414 


0.0052 


-0.0051 


(56) 
(57) 

liJ 


Radiative and 
VP corrections 
of order Q^(Za)^ 


0.0178 


0.0025 


0.0209 


0.0029 


-0.0027 


(58) 
(59) 
(63)J21] 


Recoil correction 
of order (Za)^ 


0.3009 


0.0428 


0.1781 


0.0253 


0.1053 


(64) 
\1] 


Nuclear structure 
correction of order (Za)'^ 


38.5711 


4.8214 


213.4218 


26.6825 


-152.6597 


(66) 
II.25J 


Nuclear structure 
correction of order (Za)^ 


-0.1464 


-0.0183 


-2.9384 


-0.3674 


2.4429 


(69) 
[21,25j 


Nuclear structure 
and VP correction 
of order a{Za)'^ 


0.2127 


0.0274 


1.4155 


0.1824 


-1.0478 


(70) 

(71) 


Nuclear structure 
and VP correction in 
SOPT of order a{Zay 


0.1327 


0.0135 


0.8913 


0.0898 


-0.6823 


(72) 
(73) 


Nuclear polarizability 
correction 
of order (Za)^ 


-0.1291 


-0.0161 


92.0511 


11.5064 


-80.6577 


(74) 
(75) 
[42j 
[411,42] 


HVP correction 
of order a{ZaY 


-0.0864 


-0.0108 


-0.1010 


-0.0126 


0.0128 


m 

[45] 
(77) 


Summary 
contribution 


1043927824 
~598.8893 


1043929722 
~650.1499 


1981268453 
~925.6873 


1981270452 
~980.2974 


101 003.3495 





